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ABSTRACT 

We present all X-ray and radio observations of the Type Iln supernova SN 2010jl. The X-ray observa¬ 
tions cover a period up to day 1500 with Chandra , XMM-Newton , NuSTAR , and Swift-XRT . The 
Chandra observations after 2012 June, the XMM-Newton observation in 2013 November, and most of 
the Swift-XRT observations until 2014 December are presented for the first time. All the spectra can 
be fitted by an absorbed hot thermal model except for Chandra spectra on 2011 October and 2012 
June when an additional component is needed. Although the origin of this component is uncertain, 
it is spatially coincident with the supernova and occurs when there are changes to the supernova 
spectrum in the energy range close to that of the extra component, indicating that the emission is 
related to the supernova. The X-ray light curve shows an initial plateau followed by a steep drop 
starting at day ^ 300. We attribute the drop to a decrease in the circumstellar density. The column 
density to the X-ray emission drops rapidly with time, showing that the absorption is in the vicinity 
of the supernova. We also present Very Large Array radio observations of SN 2010jl. Radio emission 
was detected from SN 2010jl from day 570 onwards. The radio light curves and spectra suggest that 
the radio luminosity was close to its maximum at the first detection. The velocity of the shocked 
ejecta derived assuming synchrotron self absorption is much less than that estimated from the optical 
and X-ray observations, suggesting that free-free absorption dominates. 

Keywords: circumstellar matter — stars: mass-loss — radiation mechanisms: non-thermal — radio 
continuum: general — Supernovae: Individual (SN 2010jl)— X-rays: general 


1. INTRODUCTION 

Type III! (narrow line) supernovae (SNe) are charac¬ 
terized by narrow emission lines atop broad wi ngs, slow 
evolu tion, and a blue continuum at early times (iSchlegell 
1199011 . Their high Ha and bolometric luminosities can 
be explained by the shock interaction of superno va (SN) 
ejecta with a dense circumstellar medium (CSM; iChugail 
1199011 . The shock waves accompanying the circumstellar 
interaction heat gas to X-ray emitting temperatures and 
accelerate particles to relativistic energies, giving rise to 
radio synchrotron emission. Indeed, Type Iln supernovae 
(SNe Iln) are among the most luminous radio and X-ray 
SNe, e.g.. SN 1986J (iBregman fc Pildislll99^. SN 1988Z 

(iFabian fc Terlevichl Il996ll. SN 1995N (iChandra et al.l 

1200511 . and SN 2006id (lCha,ndra et'aril2012bll . 

Although high CSM densities should enable radio 
and X-ray emission, few SNe Iln are detected in these 
bands. Amongst the detected ones, the X-ray and ra¬ 
dio light curves of these SNe cover a r an ge of luminosi¬ 
ties (e .g., iDwarkadas fc Gniszkol [2?)T^ . Ivan Dvk et al.l 
(|1996ll carried out a study of 10 SNe Iln with the Very 
Large Array (VLA), but did not detect radio emis¬ 
sion from any of them. Type Iln SN 1998S was not 
particularly luminous at radio and X-ray wavelengths 

^ National Centre for Radio Astrophysics, Tata Institute of 
Fundamental Research, Pune University Campus, Pune 411 007, 
INDIA, poonam@ncra.tifr.res.in 

^ Department of Astronomy, University of Virginia, P.O. Box 
400325, Charlottesville, VA 22904-4325 

® Institute of Astronomy of Russian Academy of Sciences, Py- 
atnitskaya St. 48, 109017 Moscow, Russia 

^ Oskar Klein Centre, Department of Astronomy, Stockholm 
University, AlbaNova, SE-106 91 Stockholm, Sweden 

® Smithsonian Astrophysical Observatory, 60 Garden St., MS- 
20, Cambridge, MA 02138 


(|Poolev et al.ll200^ . which can be attributed to a rel¬ 
atively low CSM density. However, SN 2006gy was very 
luminous at optical wavelengths, implying a very high 
CSM density, but was not luminous at X-ray wavelengths 
(|Ofek et al.f 120071 : iSmith et al.l 1201^ . The lack of X- 
ray emission here can be attributed to mechanisms that 
suppress the X-ray emission at high density, including 
photoelectric absorption, inverse Compton losses of hot 
shocked electrons, and Compton cooling in t he slow wind 
(jChevalier fc Irwial2012l: ISvirski et al.ll201^ . The X-ray 
luminosity of a SN may initially increase with CSM den¬ 
sity, but eventually turns over because of a variety of 
effects that suppress X-ray emission. 

In this paper, we discuss X-ray and radio observa¬ 
tions of SN 20I0jl, which may be close to the case of 
a maximum X-ray luminosity. SN 2010jl was discov¬ 
ered with a magnitude of 13.5 in unfiltered CCD im¬ 
ages with a 0.40- m reflector at Portal, AZ , U.S.A. on 
2010 November 3 (|Newton fc PuckettI 1201 011 and bright¬ 
ened to mag 12.9 over the next day, showing that it was 
discovered at an early phase. SN 2010jl is at a posi- 
tion a = 09^42( "53C337. 5 = -b09°29'42."13 (J2000) 
(jOfek et al.ll2014[j . associated with a galaxy UGC 5189A 
at a distance of 49 Mpc {z = 0.0107), implying that 
SN 2010jl belongs to the class of luminous SNe Hn with 
an absolute visual magnitude < —20. Pre-discovery 
observ ations indica,te an ex plosion date in ear ly 2010 Oc¬ 
tober (jStoll et al.l 1201 Ih . lOfek et HI (|2014[ j argue for 
an explosion date around 15 — 25 days before I-band 
maximum, i.e. around JD 2,455,469-2,455,479, or 2010 
September 29-October 9. We assume 2010 October 1 to 
be t he explosion date for SN 2010jl throughout this pa¬ 
per. iStoll et m 1)201 Ih found that the host galaxy for SN 
2010jl is of low metallicity, supporting the emerging trend 
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that luminous supernovae (SNe) occur in low metallicity 
environments. They determined the metallicity Z of the 
SN region to be < 0.3 Zq. We take this upper limit as 
the metallicity of the gas in the galaxy. 

After the X-ray Telescope (XRT) on-board Swift de- 
tected X-rays from S N 2010jl on 2010 November 5.0 —5.8 
(|Immler et al.l[20T^ . we triggered our Chandra Target of 
Opportunity (ToO) observing program in 2010 December 
an d 2011 Octobe r . Thes e observations were presented 
in iChandra et al.l (|2012af l , in which a ve ry rapid evolu¬ 
tion of the column density was reported. IChandra et al.l 
(j2012at) also reported a constant X-ray flux of the SN, 
consistent with optical w avelengths where it displayed a 
flat l i ght curve early on (Zhang et al.l l2012t lOfek et ^ 
12014 Ihdansson et al'll2014 ). While the peak R-band lu¬ 
minosity of SN 2010jl is smaller than the super-luminous 
class of supernovae, such as SN 2006gy, SN 2006tf and SN 
1997cv. SN 2010il is the most luminous X-ray SN so far. 
lOfek et ^ (12014 reported simultaneous NuSTAR and 
XMM-Newton observations and determined the temper¬ 
ature of the shock. Because of the high temperature 
of the supernova emission, the hard X-ray sensitivity of 
NuSTAR was crucial to obtain a reliable temperature es¬ 
timate. They estimated the shock velocity to be ^ 3000 
kms“^. Given the estimate of the shock velocity and the 
total luminosity of the SN, it is possible to estimate the 
density profile of the CSM if the forward shock wave is 
radiative. With this assumption the presupernova star 
lost 3 — 10 Mp in the decades prior to the explo¬ 
sion (iZhang et al.ll2012l : lOfek et ^112014 iFransson et al.l 

[Ml. 


The H lines in SN 2010jl showed a narrow component 
with an expansion velocity ~ 100 kms“^ coming from 
the CSM; along with broad wings which, at early times, 
are well fitted by an electron scattering p rofile produced 
by the thermal velocit ies of electrons (jFransson et al.l 
I20l4l^ang et aD2012[) . Here, the line profiles do not re¬ 
flect the bulk motions of the SN and the high velocity re¬ 
gions are presumably obscured by the circumstellar gas. 
However, over the first 200 days, the broad componen t 
shifts to the blue by ^ 700 kms~^ ( Fransson et al.ll2011. 

I Smith et al.l (j2012fl . iMaeda et al.l ( 20131 1. and iGall et al.l 
(|2014 1 have explained this shift to be due to the forma¬ 
tion of dust in the dense shell resulting from circumstellar 
interaction or in the freely expanding ejecta. However, 
IFransson et al.l (j20l4 1 argue that dust formation is un¬ 
likely and attribute the line shift to radiative acceleration 
of circumstellar gas; in this case, the broadening of the 
lines is due to electron scattering. Although the situa¬ 
tion with the H lines is ambiguous, there is clearer evi- 


dence for high velo city motion in the He I A10830 line. 
iBorish et al.l ([20151 1 And a blueshifted shoulder in the 
A10830 line between 100 and 200 days that is likely due to 
the ejecta emission up to a velocity of 4000—6000 kms“^, 
in rough agreement with the velocity deduced from the 
temperature of X-ray emitting gas at a later time. 

In this paper, we carry out a comprehensive analysis 
of all the X-ray and radio observations for SN 2010jl . 
The radio detection is being reported for the first time. 
In §2, we provide details of observations for SN 2010jl . 
In §3 we present analysis and interpretation of the X-ray 
emission and in §4 for the radio data. In §5 we discuss our 
main results and interpretation in view of multiwaveband 
data. The main conclusions are listed in S6. 


2. OBSERVATIONS 
2.1. X-ray Observations 

The Swift-XRT detection of SN 2010jl allowed us to 
trigger our approved Chandra Cycle 11 program and the 
first observations were made on 2010 December 7 and 
2010 December 8 for 19 and 21 ks, respectively. The ob¬ 
servations were made using the ACIS-S detector with no 
grating in the A7FAINT mode. Afterwards, we observed 
SN 2010jl on 2011 October 17 (41 ks exposure) and 2012 
June 10 (40 ks exposure) under Cycle 13 using Chandra's 
ACIS-S detector. Our most recent observation was on 
2014 June I for a 40 ks exposure under Cycle 15. We also 
observed SN 2010jl with XMM-Newton for a 52.2 ks du¬ 
ration on 2013 November 1. In addition, we use publicly 
available archival data from HEASARCQ- These were 10 
ks Chandra data on 2010 Nov 22, 12.9 ks XMM-Newton 
data observed on 2012 November 1 a nd NuSTAR dat a 
observed for 46 ks on 2013 October 5 (|Ofek et al.ll20l4) . 
as well as several Swift-XKT data sets taken between 
2010 November 5 and 2014 December 24. Table |T] gives 
details of all the X-ray observations used in this paper. 

For the Chandra data analysis, we extracted spec¬ 
tra, response and ancillary matrices using Chandra 
Interactive Analysis of Observations software (CIAO; 
iFruscione et al.l 1200611 . using task specextractor. The 
CIAO version 4.6 along with CALDB version 4.5.9 was 
used for this purpose. To extract the spectra and re¬ 
sponse matrices for XMM-Newton data, the Scientific 
Analysis System (SAS) version 12.0.1 and its standard 
commands were used. We extracted the spectra from 
NuSTAR data using The NuSTAR Data Analysis Soft¬ 
ware (NUSTARDAS) version 1.3.1. The task nupipeline 
was used to generate level 2 products and nuproducts was 
used to generate level 3 spectra and matrices. The Swift- 
XRT spectra and response matrices were e xtracted us- 
ing online XRT products building pipelin(£| ([Evans et al.l 
l2009t iGoad et M 1200711. T he HEAsofijfj package xspec 
version 12.1 ([Arnaudl 1199^ was used to carry out the 
spectral analysis. 


2.2. Radio Observations 

The radio observations of SN 2010jl were carried out 
using the Expanded Very Large Array (EVLA) tele¬ 
scope, later renamed to Karl G. Jansky Very Large Array 
(JVLA), starting from 2010 November 6 until 2013 Au¬ 
gust 10. The observations were carried out at 33 GHz 
(Ka band), 22 GHz (K band), 8.5 GHz (X band) and 
5 GHz (C band) frequency bands for 30 minute to 1 
hour durations. Each observation consisted of the flux 
calibrator 3C286 and a phase calibrator. The phase cal¬ 
ibrator was J1007-I-1356 in most cases, and J0954-I-1743 
in a few cases. The bandwidths used in the EVLA data 
and JVLA data were 256 MHz and 2048 MHz, respec¬ 
tively. In some of the EVLA observations, each 128 MHz 
subband was tuned to 4.5 and 7.5 GHz bands in order 
to estimate the flux density at the above two frequen¬ 
cies. The data were analyzed using t he Gommon As- 
tronomy Software Applications (CASA; iMcMullin et al.l 

® heasarc.gsfc.nasa.gov 

^ http://www.swift.ac.uk/user_objects/ 

® http://heasarc.gsfc.nasa.gov/docs/software/lheasoft/ 
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Table 1 

Details of X-ray observations for SN 2010jl 


Date of 
Observation 

Mission 

Instrument 

Obs. 

ID 

Exposure 

ks 

^^exp 

days ^ 

2010 Nov 05.02 

Swift 

XRT 

00031858001 

1.38 

36.0 

2010 Nov 05.08 

Swift 

XRT 

00031858002 

1.98 

36.1 

2010 Nov 05.67 

Swift 

XRT 

00031858003 

6.96 

36.7 

2010 Nov 05.88 

Swift 

XRT 

00031858004 

4.81 

36.9 

2010 Nov 06.08 

Swift 

XRT 

00031858005 

1.97 

37.1 

2010 Nov 07.02 

Swift 

XRT 

00031858006 

2.33 

38.0 

2010 Nov 08.16 

Swift 

XRT 

00031858007 

2.80 

39.2 

2010 Nov 09.02 

Swift 

XRT 

00031858010 

1.68 

40.0 

2010 Nov 09.10 

Swift 

XRT 

00031858008 

0.18 

40.1 

2010 Nov 09.10 

Swift 

XRT 

00031858009 

0.58 

40.1 

2010 Nov 11.06 

Swift 

XRT 

00031858011 

2.30 

42.q 

2010 Nov 12.03 

Swift 

XRT 

00031858012 

0.56 

43.0 

2010 Nov 12.03 

Swift 

XRT 

00031858014 

11.13 

43.0 

2010 Nov 12.16 

Swift 

XRT 

00031858013 

2.09 

43.2 

2010 Nov 13.78 

Swift 

XRT 

00031858015 

2.22 

44.8 

2010 Nov 14.11 

Swift 

XRT 

00031858016 

2.41 

45.1 

2010 Nov 15.58 

Swift 

XRT 

00031858017 

1.83 

46.6 

2010 Nov 16.71 

Swift 

XRT 

00031858018 

2.29 

47.7 

2010 Nov 17.45 

Swift 

XRT 

00031858019 

2.13 

48.5 

2010 Nov 20.07 

Swift 

XRT 

00031858020 

2.11 

51.1 

2010 Nov 22.03 

Chandra 

ACIS-S 

11237 

10.05 

53.0 

2010 Nov 23.01 

Swift 

XRT 

00031858021 

2.41 

54.0 

2010 Nov 26.16 

Swift 

XRT 

00031858022 

2.41 

57.2 

2010 Nov 29.69 

Swift 

XRT 

00031858023 

2.27 

60.7 

2010 Dec 02.13 

Swift 

XRT 

00031858024 

2.09 

63.1 

2010 Dec 05.67 

Swift 

XRT 

00031858025 

2.42 

66.7 

2010 Dec 07.18 

Chandra 

ACIS-S 

11122 

19.05 

68.2 

2010 Dec 08.03 

Chandra 

ACIS-S 

13199 

21.05 

69.0 

2011 Apr 24.56 

Swift 

XRT 

00031858026 

7.89 

206.6 

2011 Apr 28.04 

Swift 

XRT 

00031858027 

2.28 

210.0 

2011 Oct 17.85 

Chandra 

ACIS-S 

13781 

41.04 

382.9 

2012 Jun 10.67 

Chandra 

ACIS-S 

13782 

40.07 

619.7 

2012 Oct 05.98 

NuSTAR 

FPMA 

40002092001 

46.11 

737.0 

2012 Oct 05.98 

NuSTAR 

FPMB 

40002092001 

46.07 

737.0 

2012 Oct 07.02 

Swift 

XRT 

00080420001 

2.65 

738.0 

2012 Oct 21.31 

Swift 

XRT 

00032585001 

8.08 

752.3 

2012 Nov 01.63 

XMM-Newton 

EPIC-PN 

0700381901 

4.04 

763.6 

2012 Nov 01.63 

XMM-Newton 

EPIC-MOSl 

0700381901 

10.11 

763.6 

2012 Nov 01.63 

XMM-Newton 

EPIC-MOS2 

0700381901 

9.73 

763.6 

2013 Jan 21.10 

Swift 

XRT 

00032585002 

8.19 

844.1 

2013 Feb 10.00 

Swift 

XRT 

00032585003 

4.78 

864.0 

2013 Feb 20.61 

Swift 

XRT 

00032585004 

5.90 

874.6 

2013 Mar 04.49 

Swift 

XRT 

00046690001 

0.92 

886.5 

2013 Mar 29.27 

Swift 

XRT 

00032585005 

18.63 

911.3 

2013 May 14.68 

Swift 

XRT 

00032585006 

6.89 

957.7 

2013 May 15.34 

Swift 

XRT 

00032585007 

5.47 

958.3 

2013 May 19.95 

Swift 

XRT 

00032585008 

4.82 

963.0 

2013 May 21.48 

Swift 

XRT 

00032585009 

4.76 

964.5 

2013 Jun 28.00 

Swift 

XRT 

00032585010 

8.20 

1002.0 

2013 Jun 28.34 

Swift 

XRT 

00032585011 

6.29 

1002.3 

2013 Nov 01.67 

XMM-Newton 

EPIC-PN 

0724030101 

52.30 

1128.7 

2013 Nov 01.67 

XMM-Newton 

EPIC-MOSl 

0724030101 

52.30 

1128.7 

2013 Nov 01.67 

XMM-Newton 

EPIC-MOS2 

0724030101 

52.30 

1128.7 

2013 Dec 11.54 

Swift 

XRT 

00032585012 

2.07 

1168.5 

2013 Dec 18.00 

Swift 

XRT 

00031858013 

0.95 

1175.0 

2013 Dec 19.07 

Swift 

XRT 

00032585014 

1.21 

1176.1 

2013 Dec 20.00 

Swift 

XRT 

00032585015 

3.20 

1177.0 

2013 Dec 24.40 

Swift 

XRT 

00032585016 

6.68 

1181.4 

2013 Dec 30.67 

Swift 

XRT 

00032585017 

0.59 

1187.7 

2014 May 11.12 

Swift 

XRT 

00032585018 

5.07 

1319.1 

2014 May 13.04 

Swift 

XRT 

00032585019 

1.02 

1321.0 

2014 May 14.64 

Swift 

XRT 

00032585020 

2.24 

1322.6 

2014 May 16.84 

Swift 

XRT 

00032585021 

3.10 

1324.8 

2014 May 18.78 

Swift 

XRT 

00032585022 

0.58 

1326.8 

2014 May 20.44 

Swift 

XRT 

00032585023 

0.32 

1328.4 

2014 Jun 01.25 

Chandra 

ACIS-S 

15869 

40.06 

1340.3 

2014 Jun 27.90 

Swift 

XRT 

00046690002 

0.92 

1366.9 

2014 Nov 30.76 

Swift 

XRT 

00032585024 

1.27 

1522.7 

2014 Dec 04.22 

Swift 

XRT 

00032585026 

1.04 

1526.2 

2014 Dec 05.08 

Swift 

XRT 

00032585027 

1.51 

1527.1 

2014 Dec 08.54 

Swift 

XRT 

00032585028 

3.35 

1530.5 

2014 Dec 09.01 

Swift 

XRT 

00032585029 

2.39 

1531.0 

2014 Dec 10.00 

Swift 

XRT 

00032585030 

0.35 

1532.0 

2014 Dec 18.51 

Swift 

XRT 

00032585031 

2.84 

1540.5 

2014 Dec 19.52 

Swift 

XRT 

00032585032 

1.52 

1541.5 

2014 Dec 24.50 

Swift 

XRT 

00032585033 

2.92 

1548.5 


^ Assuming 2010 October 1 to be SN 2010jl explosion date. 
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I2007f) . The VLA Calibration pipeline^ was used for flag¬ 
ging and calibration purposes. However, in several cases, 
extra flagging was needed. In those cases, flagging and 
calibration was done manually. The images were made 
with CASA task ‘clean’ in which “briggs” weighting with 
a robustness parameter of 0.5 was used. For the 2 GB 
bandwidth data, “mfs” spectral gridding mode with two 
Taylor coefficients was used to model the sky frequency 
dependence. The observational details are presented in 
Table [H 

3. X-RAY analysis AND INTERPRETATION 

3.1. Analysis of the Contaminating Sources 

SN 2010jl host galaxy belongs to UGC 5189 group 
of galaxies, which has a size of 1.7' centered at a = 
09''42™54."72, <5 = -P09°29'01."4 (J2000). A NASA Ex- 
tragalactic Database^ (NED) search shows that there 
are three sources within 3" of the SN 2010jl posi¬ 
tion. These are UGG 05189 NEDOl (UGG 5189A, 
a = 09''42™53A434, <5 = -K09°29'41."87 (J2000)), 
MCG -P02-25-021 GROUP (a = 09''42™53'’, ^ = 
-f09°29.'7 (J2000)) and SDSS J094253.47-P092943.5 
{a = 09''42™53*.47, S = -P09°29'43."51 (J2000)) at dis¬ 
tances of 1.44", 2.34" and 2.46" away from the SN, re¬ 
spectively. The SDSS J094253.47-P092943.5 and UGG 
05189A sources have been identified as galaxies, whereas 
MGG -1-02-25-021 GROUP is a group of galaxies. Since 
there are no X-ray archival data at the SN 2010jl field of 
view (FoV), we could not ascertain whether these three 
sources were X-ray emitters or not. For this reason we 
started our analysis with the Chandra data. Since Chan¬ 
dra has excellent spatial resolution, and can separate out 
the nearby sources. 

The SN 2010jl FoV in Chandra observations at various 
epochs show UGC 5189A to be an X-ray emitter. How¬ 
ever, no X-ray emission in seen from MCG -1-02-25-021 
GROUP or SDSS J094253.47-P092943.5. Thus we need 
to make sure that the UGG 5189A does not contaminate 
the SN 2010jl flux in Chandra data. We define three 
boxes in the Chandra FoV as shown in the left panel of 
Fig. [B Box A of size 2.3" x 3.2" covers SN 2010jl, while 
box B of size 2.0" x 2.5" covers UGG 5189A. We also 
extract a 4.3" x 4.0" box covering both SN 2010jl and 
UGG 5189A centered at (J2000) a = 09''42™53."361 
S = -|-09°29'41."55 (Box C). The background region is 
chosen in a source free area with a 9.0" x 9.0" box. 

In order to estimate the contamination from UGC 
5189A , we start our analysis with the Chandra observa¬ 
tions on 2010 December 7 and 8. We extract the spectra 
of the SN alone from Box A, and of UGC 5189A from 
box B. We also extract the combined spectrum from Box 
C. The spectra are grouped into 15 channels for Boxes A 
and C and x^-statistics is used. However, due to a small 
number of counts in the Box B, the spectrum was binned 
into 5 channels per bin and C-statistics were used to fit 
the data. We fit the UGC 5189A spectrum with an ab¬ 
sorbed power law model, whereas, we use the Astroph ys- 
ical Plasma Emission Code (apec; ISmith et al.ll200i1l to 
fit the spectrum for SN 2010jl. The apec gives a fit to an 
emission spectrum from collisionally-ionized diffuse gas. 


The parameters of this model are the plasma tempera¬ 
ture, metal abundances and redshift. In the spectrum 
from Box C, we ht the joint absorbed thermal plasma 
and absorbed power law spectra. Here we fix the re¬ 
spective parameters to the best fit values obtained from 
the spectral fits of Box A and Box B; however, we let 
the normalizations vary. We estimate the 0.2-10 keV 
fluxes in Boxes A, B, and C to be (6.58 ± 0.38) x 10“^^ 
erg cm“^ s“^, (1.50 ± 0.23) x 10“^^ erg cm“^ s“^ and 
(6.93 ± 0.32) X 10“^^ erg cm“^ s“^, respectively. The 
flux in Box C matches the total flux from Boxes A and B 
within the error bars. In addition, the UGC 5189A flux 
is 40 times weaker than the SN flux. Thus the contami¬ 
nation to the SN flux due to UGC 5189A is insignificant 
in the Chandra data. 

Now, we carry out a joint fit to the spectra from UGC 
5189A (Box B) at all five Chandra epochs of observation 
with an absorbed power law model. We assume that the 
absorption column density and power law index do not 
change at various epochs and enforce these parameters 
to be the same at all epochs by linking them in the fits. 
However, we let normalizations vary independently to ac¬ 
count for variable X-ray emission. Since the counts are 
few, we use C-statistics to fit the data. The model is 
best fitted by an absorbed power law with an absorption 
column density of Nh = ( 1 . 82 ti; 63 ) x 10^^ cm“^ and a 
photon power law index F = 1 . 161 q; 4 q. This type of pho¬ 
ton index is consistent if the X-ray emission in the galaxy 
is mainly from a collection of X-ray binaries (XRBs) or 
from a combination of XRBs and diffuse gas. The best 
fit C-statistic is 22.94 for 30 degrees of freedom. The ab¬ 
sorption column density is much higher than the Galac¬ 
tic absorption which is A^j^ (Galactic) = 3 x 10^° cm~^. 
The remaining column density A^//(Host) = 1.52 x 10^^ 
cm“^ must be coming from the host galaxy UGC 5189A. 
To take care of the low metallicity of the host galaxy, 
we refit the data with an absorbed power law with col¬ 
umn density Nh to be Nh = Wjj (Galactic) -I-A^^r(Host). 
We use solar metallicity for Wjj (Galactic) and fix it to 
3 X 10^° cm“^. We fix the metallicity of Nh (Rost) to 
be 0.3 solar and treat A^//(Host) to be a free parameter. 
The best fit values are Nh (Rost) = (d.lOtl;^) ^ 

cm“^ and F = l.ldj'jQ'gg. The increase in Vi/(Host) 
is due to the lower value of metallicity for the host 
galaxy. This simply means that the equivalent hydro¬ 
gen column density has to be 1/0.3 times larger to ac¬ 
count for the same absorption of X-rays by metals in a 
0.3 Zq metallicity environment. Our attempt to fit the 
UGC 5189A spectra by fixing the column density to that 
of the Galactic value results in a relatively flat power 
law index F = 0.85lg'2Q, which is nonphysical. There 
is an additional evidence of higher neutral HI column 
density towards UGC 5189 from the Giant Metrewave 
Radio Telescope (GMRT) 21 cm radio data obtained in 
2013 Nov-Dec (|Chengalur et al.l 1201411 . At a position 
a = 09''42"'53".434, S = -L09°29'41."87 (J2000), they 
find the HI flux to be 5.8 mJy kms“^, which translates 
to a HI column density of 2.4 x 10^^ cm“^. 

We list the absorbed and unabsorbed fluxes and un¬ 


absorbed luminosities of UGC 5189A in the 0.2-10 keV 
https://science.nrao.edu/facilities/vla/data-processing/pipefHy&® five Chandra epochs in Table[2] The flux varies 
J http://ned.ipac.caitech.edu/ at most by a factor of 1.8 in these epochs, so we re- 
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Table 2 



0.2- 

TO keV X-ray fluxes of UGC 5189A in Chandra 

observations 


Date of 
Observation 

^I'exp 

days^ 

Count 

Rate 

Abs. Flux 
erg cm“^ s“^ 

Unabs. Flux 
erg cm“^ s~^ 

Unabs. Luminosity 
erg s“^ 

2010 Nov 22.03 

2010 Dec 07.18-8.03 

2011 Oct 17.85 

2012 Jun 10.67 

2014 Jun 01.25 

Joint fit 

53 

68.2-69.0 

382.9 

619.7 

1340.3 

(1.50 ±0.39) X 10-3 
(0.87 ±0.15) X 10-3 
(1.07 ±0.17) X 10-3 
(1.37 ±0.19) X 10-3 
(0.92 ±0.15) X 10-3 

(2.93 ±0.66) X lO-i'^ 
(1.61 ±0.25) X lO-i'^ 
(1.99 ±0.28) X lO-i'^ 
(2.63 ±0.33) X 10-14 
(1.69 ±0.26) X 10-14 
(2.02 ±0.14) X 10-14 

(3.23 ± 0.73) X 10-14 
(1.77 ±0.28) X 10-14 
(2.19 ± 0.31) X 10-14 
(2.90 ± 0.37) X 10-14 
(1.86 ± 0.28) X 10-14 
(2.22 ± 0.16) X 10-14 

(9.27 ±2.09) X 1039 
(5.10 ±0.79) X 1039 
(6.29 ±0.90) X 1039 
(8.34 ± 1.05) X 1039 
(5.34 ±0.81) X 1039 
(6.13 ±0.44) X 1039 


^Jot©. — T1 i 6 fluxGS cir© d©riv©d cis d©t3.il©d in El] 

^ Assuming 2010 Octob©r 1 to b© SN 2010jl ©xplosion dat©. 


Table 3 

0.2-10 keV fluxes of 6 sources (Fig. ^ right panel) within 21" radius of the SN 2010jl position in Chandra observations 


Date of 
Observation 

^texp 

days^ 

Count 

Rate 

Abs. Flux 
erg cm“^ s~^ 

Unabs. Flux 
erg cm~^ s~^ 

Unabs. Luminosity 
erg s“^ 

2010 Nov 22.03 

2010 Dec 07.18-8.03 

2011 Oct 17.85 

2012 Jun 10.67 

2014 Jun 01.25 

Joint fit 

53 

68.2-69.0 

382.9 

619.7 

1340.3 

(2.28 ±0.58) X 10-3 
(1.68 ±0.26) X 10-3 
(3.53 ±0.34) X 10-3 
(2.23 ±0.30) X 10-3 
(2.82 ±0.28) X 10-3 

(2.17 ±0.45) X 10-44 
(1.51 ±0.19) X 10-14 
(3.17 ±0.27) X 10-14 
(2.27 ±0.24) X 10-14 
(2.56 ±0.24) X 10-14 
(2.22 ±0.11) X 10-14 

(3.19 ± 0.66) X 10-14 
(2.21 ± 0.28) X 10-14 
(4.64 ±0.40) X 10-14 
(3.33 ±0.35) X 10-14 
(3.75 ±0.36) X 10-14 
(3.30 ±0.17) X 10-14 

(8.81 ± 1.82) X 1039 
(6.10 ±0.76) X 1039 
(12.84 ±1.10) X 1039 
(9.21 ±0.95) X 1039 
(10.37 ±0.99) X 1039 
(9.11 ±0.47) X 1039 

Note. — The fluxes 

are derived 

as detailed in ^3.11 





^ Assuming 2010 October 1 to b© SN 2010jl ©xplosion dat©. 



53.6 53.5 53.4 9 : 42 : 53.3 53.2 53.1 

Right ascension 



55.0 54.5 9:42:54.0 53.5 53.0 52.5 52.0 


Right ascension 


Figure 1. Left Panel: SN 2010jl field of view (FoV) from Chandra observations. The left green box (2.0" x 2.5") marked ‘B’ is used 
to extract the spectrum for UGC 5189A. The right green box (2.3" x 3.2") marked ‘A’ has been used to extract the SN 2010jl spectrum. 
The largest red box (4.3" x 4.0") marked ‘C’ includes both SN 2010jl and the UGC 5189A. Right Panel: Field of view from Chandra 
observations. Here SN 2010jl is the brightest source at the top part of the image center. The green circles are the 6 nearby sources within 
a 21" radius centered at the SN 2010jl position, excluding UGC 5189A. 
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fit the data assuming constant flux at all the epochs. 
This also gives a reasonable fit with a best fit C-statistic 
of 31.4 for 36 degrees of freedom. Here the best fit 
values are iV//(Host) = x 10^^ cm“^ and 

r = l.lSto^sg. The 0.2-10 keV absorbed (unabsorbed) 
flux of UGC 5189A is (2.02 ± 0.14) x 10“^^ erg cm“^ s“^ 
((2.22 ± 0.16) X 10“^'^ erg cm“^ s“^), which translates 
to an unabsorbed luminosity of (6.13 ± 0.44) x 10^® erg 
s“^ (Tabled). We also plot the contour levels of the best 
fit column density and photon index in Figure [2l We 
note that there is a large uncertainty in the absorption 
column density. 

The spatial resolution of the Swift-XRT , XMM- 
Newton and NuSTAR observations are not as good as 
that of Chandra , so we have to take care of contami¬ 
nation from more distant sources while analyzing these 
data. We looked for the contaminating sources within 
a 60" region centered at SN 2010jl. As shown in the 
Chandra FOV (right panel of Fig. [1]), there are 6 sources 
within 21" of the SN position, in addition to UGC 5189A. 
There are no additional sources between 21" and 60" ra¬ 
dius centered at the SN. We extracted spectra in the 
0.2-10 keV range from Chandra data at each epoch 
for these 6 sources and carried out a joint fit, assum¬ 
ing the flux did not change at various Chandra epochs. 
Their spectra are best fit with a column density of 
NniHost) = (4.72l3;gg) x 10^^ cm“^ and a power law in¬ 
dex of F = 2.051 q 37 (reduced = 1.07). We note that 
the absorption column density is similar to that obtained 
for UGC 5189A, re-confirming that the host galaxy UGC 
5189 contributes a significant amount of X-ray absorbing 
column. The 0.2-10 keV absorbed (unabsorbed) flux of 
these sources combined is (2.22 ±0.11) x 10“^^ erg cm“^ 
s“^ ((3.30 ± 0.17) X 10“^"^ erg cm“^ s“^), which trans¬ 
lates to an unabsorbed luminosity of (9.11 ±0.47) x 10^® 
erg s“^. We also attempted to carry out fits where we 
let normalizations at each epoch vary. The absorbed flux 
changed at most by a factor of 2, the reduced =0.84 
improved significantly. In Table[3]we give the 0.2—10 keV 
fluxes of these sources. 

3.2. Analysis of SN 201 Ojl 

iChandra et al.l (j2012al l analyzed Chandra data from 
2010 December and 2011 October observations. In both 
cases, their best fit temperature values always hit the 
hard upper limit of the models in XSPEC. The power 
law models were disca rded since they gave a n unphys- 
ically hard spectrum. IChandra et al.l (I2012af) preferred 
a thermal model, noting that the plasma giving rise to 
the X-ray emission is sufficiently hot that Chandra is not 
sensitive to the high plasma temperature. However, we 
now have the advantage of having SN 2010jl observa¬ 
tions with NuSTAR which has sensitivity in the range 
3-80 keV. Because of the above complications, we use 
NuSTAR data to determine the shock temperature and 
then use the same temperature for the analysis of the 
rest of the data. 

3.2.1. NuSTAR data 

NuSTAR observed SN 2010jl on 2012 October 5 
(lOfek et all [Ml . SN 2010jl was also observed with 
XMM-Newton on 2012 November 1. To get spectral cov- 


68%, 90% and 99% Confidence contours for Photon index and 



Figure 2. The 68% (black), 90% (red) and 99% (green) confidence 
contours for best fit column density {N^) and power law photon 
index (F) for the UGC 5189A X-ray spectra obta ined from the joint 
fit of Chandra data at various epochs (see gn}. The uncertainty 
in the column density is large. 

erage in a wider energy range, we carried ont a joint 
fit to both NuSTAR and XMM-Newton spectra. We 
used the best fit parameters of the 2012 June 1 Chan¬ 
dra data for UGC 5189A and the other 6 contaminat¬ 
ing sources to subtract out their contamination to SN 
flux (see 113.11) . For the XMM-Newton observations, we 
used spectra from MOSl, MOS2 and PN CCD arrays, 
whereas for the NuSTAR observations, we used spectra 
from both FPMA as well as FPMB focal plane modules. 
The total absorption column density in the models are 
Nh, where Nh = V//(Galactic) ±Nh(H ost) ±Nh(CSM). 
Here V//(Galactic) and Nh(H ost) are explained in 113.11 
and Nh {GSM) is the column density due to the SN 
CSM. For V^f(CSM), we fix the metallicity to be 0.3 
solar. The data are best fit with an absorbed thermal 
plasma model. The reduced x^ is 0.97 for 165 degrees 
of freedom. The best fit column density for SN 2010jl is 
Nh{CSM) = {G.GTfi'g'^) X 10^^ cm“^ and the plasma 
temperature is kT = IS.OO^lge keV. This temperature 
is consi stent with that f ound b v lOfek et al.l l)2ni4[ ). 

Since IChandra et al.l (|2012afl claimed the presence of 
a 6.33 keV Fe Ka line in their December 2010 Chan¬ 
dra spectrum, we also attempted to add a Gaussian com¬ 
ponent around the same energy and refit the spectra. 
There was no significant change in the quality of the fit. 
Thus Fe K-a line is not significant here. In Figured we 
plot the best fit model as well as the contour diagram of 
Nh {CSM) versus kT. The 0.2-80 keV absorbed (unab¬ 
sorbed) flux of the SN is (5.24 ±0.25) x 10“^^ erg cm“^ 
s-i ((6.01 ±0.28) X 10-13 erg cm'^ s"!). The 0.2-10 keV 
absorbed (unabsorbed) flux of the SN is listed in Table 

El 
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Nov 2012 XMM spectra of SN 201 Ojl 
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68%, 90% and 99% contours for N^ & kT for NuSTAR+XMM data 



68%, 90% and 99% contours for N^ for Nov 2012XMM data 



0.6 

Nh(10^^ 


0.8 

cm-2) 


68%, 90% and 99% contour plots for N^ for Nov 2013 XMM data 



0.25 

N„ (10^= 


Figure 3. The spectra of SN 2010jl and their best fit models as described in § im] Upper panel: The best fit model for the 2012 October 6— 
November 1 joint XMM-Newton (green: PN, black: MOSl, red: MOS2) and NuSTAR (blue: FPMA, cyan: FPMB) spectra. The data are best 
fit with an absorbed apec model. The right plot shows the 68% (black), 90% (red) and 99% (green) confidence contours for the best fit CSM 
column density {Nh) and the plasma temperature {kT). Middle panel: The spectrum and column density confidence contours for 2012 Nov 1.63 
XMM-Newton data. Lower panel: The spectrum and column density confidence contours for 2013 Nov 1.67 XMM-Newton data. In all the cases, 
Nfj is well constrained. The colors in the middle and lower panels have the same association as explained in the top panel. Here the normalized 
counts s~^ kev~^ in the y-axes in the left plots of the spectra are normalized counts s“^ kev~^. The word normalized indicates that this plot has 
been divided by the effective area, the value of the EFFAREA keyword, in the response file associated with each spectrum. In all the left side plots, 
the residuals are in terms of as with error bars of size one. 
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Table 4 

Best fit models at various epochs 


Date of 

Obsn 

Instrument 

Model 

Param-1 

Nh (cm-2) 

Param-2 

Nh3 (cm-2) 

Param-3 

Reduced 

2010 Nov 5-20 

Swift-XRT 

E Nh* * Ti + * {kT + Gauss) 

i=l 

(17.67l7®gf ) ^ 


Gauss= 6.391 q23 

1.730 

2010 Nov 22.03 

Chandra 

Nfj * kT 

(9.59+^;';^) X 10^“ 



1.570 

2010 Nov 23-Dec 05 

Swift-XRT 

E Nh* * r, + * fcT 

i=l 

9.59 X 1023 (fixed) 



1.350 

2010 Dec 7-8 

Chandra 

N^ * {kT + Gauss) 

(9.47t";“) X 1023 


Gauss= 6.33llog 

1.90 

2011 Apr 24-28 

Swift-XRT 

E Nh* *ri + Nfj*kT 

i = l 

i5.7i+ill) X 10^3 



O. 92 O 

2011 Oct 17.85 

2011 Oct 17.85 

2011 Oct 17.85 

Chandra 

Chandra 

Chandra 

Nfj * kT 

N^ * kT + Nh^ * r 

N^ * kT + Nh^ * r 

""(1.63) X 1023 
(2.05j:°-24) X 

(2.89t°J6) X ^(,23 

4 X 102l (fixed) 
(3.00llJ^) X 1021 

r = 1.7 (fixed) 
r = 1.7 (fixed) 

2.07 

1.37 

1.11 

2012 Jun 10.67 

2012 Jun 10.67 

2012 Jun 10.67 

Chandra 

Chandra 

Chandra 

N^A^kT 

* kT + Nh 3 * r2 

N^ * kT + Nh 3 * r2 

"‘(3.74 X 1022) 

(l.llt°-2|) X 1023 

(1.09t°;23) X 1023 

0.4 X 1021 (fixed) 
(3.70im) X 1021 

r 2 = 1.7(fixed) 
r 2 = 1.7 (fixed) 

2.71 

0.93 

0.94 

2012 Oct 7-21 

Swift-XRT 

E Nh* *ri + Nfj*kT 

i = l 

6.67X1021 (fixed) 



2.860 

2012 Oct 5-Nov 1 

NuSTAR ,XMM-Newton 

E Nh* *ri + Nfj*kT 

i = l 

(6.6712-47) X 1021 



0.97 

2013 Jan 5-Mar 29 

Swift-XRT 

E Nh* *ri + Nfj*kT 

i = l 

(4.0313-32) X 1021 



1.480 

2013 May 14-Jun 28 

Swift-XRT 

E Nh* *ri + Nfj*kT 

i = l 

(5.581111) X 10^^ 



I. 87 O 

2013 Nov 1.67 

XMM-Newton 

E Nh* *ri + Nfj*kT 

i = l 

(2.64l° ®®) X 1021 



1.41 

2013 Dec 11-30 

Swift-XRT 

E Nh* *ri + N*^*kT 

i = l 

2.64 X 102l (fixed) 



I. 22 O 

2014 Jun 1.25 

Chandra 

Nfj * kT 

(6.8211"^) X 1021 



1.24 

2014 May 11-Jun 27 

Swift-XRT 

E Nh* *ri + Nfj*kT 

i = l 

6.82 X 1021 (fixed) 



O. 65 O 

2014 Nov 30-Dec 24 

Swift-XRT 

E Nh* *ri + Nfj*kT 

i=l 

(1.41112^90) X 1021 



O. 43 O 


Note. — Except for the joint XMM-Newton and NuSTAR spectrum, which gave us the best fit temperature of kT = IS.QQl^'gg keV, everywhere 
else the temperature has been kept fixed to a value of 19 keV. The kT corresponding to all the models in the model corresponds to an apec thermal 
plasma model. The i = 1 index in Nni is for UGC 5189A and i = 2 for 6 nearby sources within a 21" error circle centered at the SN position. The 
parameters with no suffix correspond to the main SN component and parameters with sufhx 3 correspond to the extra soft component present in the 
2011 October and 2012 June data. We have used Nh* ~ S x 10^° + Nni, Nh* = 3 x 10^° + Nh and Nhs = 3 x 10^° + Nh 3 - This is to account for 
the contribution from Galactic absorption whose metallicity is fixed to solar. For host and GSM contributions, the metallicity is fixed to 0.3 solar. For 
data with multiple fits, the models in bold are considered to be the models best representing the respective spectra. 

^ Since > 2, XSPEG did not calculate the error. 

^ G-statistics were used to fit the data. Equivalent Sire mentioned. 
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November 2010 Chandra spectrum of SN 201 Oil 



Figure 4. The best fit models to Chandra spectra of SN 2010jl at 
2010 Nov, 2010 Dec and 201 4 June, respectively. The best fit 
models are explained in § 13.21 and listed in Table U 


3.2.2. Chandra data 

For all data, we fix the temperature to be kT = 19 keV 
(see il3.2.1l) . Even though this temperature is probably a 
lower limit to the temperature for the SN shock at ear¬ 
lier epochs, it is the best available temperature estimate. 
This will probably introduce some errors in the SN flux 
and the column density estimates. However, we discuss 
in aSD that the uncertainties due to the assumption of 
constant temperature are not signihcantly large. 

The Chandra data obtained in 2010 November 22 have 
very few counts, so we group the spectrum into 5 counts 
per bin and use C-statistics to fit the data. We fix all 
the parameters except the normalization and Nh- For 
a metallicity of 0.3, the data are best fit with a column 
density of iViy(CSM) = (9.59t|®°) x 10^3 cm-^. 

The 2010 December 7-8 spectra were grouped in 15 
counts per bin and y^-statistics were applied to obtain 


best fits. The spectra are best fit with a column density 
of Nh{CSM) = (9.47to;52) X 10^^ cm“^. However, there 
is an indication of extra emission com ponent around 
6 keV energy, which we also seen by IChandra et al.l 
(|2012ali and was associated with Fe K-a line. In our 
current fits, the line is best fit with a Gaussian of en¬ 
ergy FiGauss = O.SS^o^os keV and width 0.19 keV. The 
0.2-10 keV unabsorbed flux in this line component is 
(4.59 ± 2.53) X 10-34 erg cir-^ s-^. 

The analysis of 2011 October data show that the 6.33 
keV iron line is not present in the spectrum, and the 
CSM column density has now decreased to a value of 
Nh{CSM) = 1.63 X 10^3 cm”^. A decrease of the col¬ 
umn density with time is expected as the shock moves 
to larger radii. However, the apec model does not ht the 
SN spectrum well (reduced = 2.07 for 73 degrees of 
freedom, also see Fig. 0). There appears to be an ex¬ 
tra component at the lower energy end of the spectrum. 
We explore three possibilities for this component. First, 
the component may be coming from the same region as 
the harder X-ray emission, which is most probably the 
forward shock. In this case, the column density for the 
soft extra component, A^//(Soft), should be the same as 
that of the 19 keV component. The second possibility 
is that this component is arising from the reverse shock. 
In this case the absorption for this component should be 
higher than that of the 19 keV component as the cool 
shell will contribute to an additional absorption. In the 
third case, we let the Nh vary independently. The first 
and second possibilities seem unlikely because by fixing 
the column density to either that of the 19 keV com¬ 
ponent, or declaring the column density associated with 
the 19 keV component to be the lower limit for this extra 
component (cool shell origin of the component), neither 
the apec nor the power law models give a good fit. The 
best fit models result in an extremely low temperature 
<0.1 keV or negative power law index, and give a 7 or¬ 
ders of magnitude higher normalization than the harder 
component, which is nonphysical. When we fit the spec¬ 
tra with a power law index of F = 1.7 and fix the column 
density for this extra component to be that of the host 
galaxy, i.e. 4.1 x lO^^ cm”^, the reduced y^ improved 
from 2.07 to 1.37. Allowing TV//(Soft) to vary freely gives 
a best fit with F7//(Soft) = 3.00l)'2o x 10^^ cm”^ and the 
reduced y^ improves significantly to 1.11. In the case 
of fixing the column density to the host galaxy absorp¬ 
tion, the 0.2-10 keV absorbed (unabsorbed) flux of the 
extra component is (3.30 ± 0.67) x lO-^^ erg cm”^ s-^ 
((4.14 ± 1.84) X 10-34 erg cm-^ §-3). In the case where 
we let the column density be a free parameter, the 0.2-10 
keV absorbed (unabsorbed) flux of the extra component 
is (1.22±0.16) X 10-33 erg cm-^ s-3 ((2.01±1.27) x 10-33 
erg cm-^ s~^)- 

Now we carry out an analysis of the 2012 June Chan¬ 
dra spectrum. Here the column density is best fit with 
a value of Nh(CSM) = 3.74 x 10^^ cm-^. The low tem¬ 
perature feature seen in the 2011 October data is promi¬ 
nent here as well and the reduced y^ for the absorbed 
thermal plasma is quite large (y^ = 2.71 for 80 degrees 
of freedom). For this additional low temperature com¬ 
ponent, we explore the same three possibilities as dis¬ 
cussed in the above paragraph for the case of 2011 Oc¬ 
tober data. In the first case, when A^^f(Soft) is fixed to 
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June 2012 Chandra spectrum of SN 201 Ojl 



Figure 5. October 2011 (left panels)and June 2012 (right panels) Chandra spectra of SN 2010jl. These are the two epochs where the evidence 
of an extra component is present at the lower energy end of the spectra. Left Panels: The Chandra best fit spectra for the 2011 October data. 
The top row is fit with only a thermal plasma model (black dashed line). The middle row includes an additional absorbed power law model (blue 
dashed line) and fixing the absorption of the extra component to that of host galaxy’s absorption, i.e. 4 X 10^^ cm~^. The lower row is the same as 
the middle row, except the column density has been kept as a free parameter. Right Panels: Same as left, but for SN 2010jl spectra of 2012 June 
Chandra data. 
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November 2010 Swift spectrum of SN 201 Ojl 


April 2011 Swift spectrum of SN 2010jl 
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May-June 2013 Swift-XRT spectrum of SN 201 Ojl 







Figure 6 . The best fit Swift-XRT spectra of SN 2010jl at various epochs. The best fit models plotted here are detailed in ^ I3.2l and listed 
Table [4| The black dashed line is the best fit apec thermal plasma model. The red line is the contribution of UGC 5189A and the blue 
line is the contribution from six nearby sources. 
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Table 5 

SN 2010jl 0.2—10 keV flux at various epochs 


Date of 
Observation 

Telescope 

^^exp 

days 

Abs. Flux 
erg cm“^ s“^ 

Unabs. Flux 
erg cm“^ s“^ 

Unabs. Luminosity 
erg s~^ 

2010 Nov 5-20 

Swift-XRT 

43.55 ± 7.53 

(4.67 ±0.81) X 10-13 

(29.54 ±7.95) x 10-13 

(8.49 ± 2.29) X lO-ii 

2010 Nov 22.03 

Chandra 

53.03 

(7.07 ± 1.15) X 10-13 

(30.56 ±4.95) X 10-13 

(8.78 ± 1.42) X lO-ii 

2010 Nov 23-Dec 5 

Swift-XRT 

60.34 ±6.33 

(5.42 ± 1.81) X 10-13 

(24.33 ±8.13) X 10-13 

(6.99 ± 2.33) X 10^1 

2010 Dec 7-8 

Chandra 

68.61 ±0.43 

(7.03 ±0.52) X 10-13 

(29.40 ±2.48) X 10-13 

(8.45 ± 0.71) X 10^1 

2011 Apr 24-28 

Swift-XRT 

208.30 ± 1.74 

(10.70 ±2.39) X 10-13 

(34.68 ±7.76) x 10-13 

(9.96 ± 2.23) X 10^1 

2011 Oct 17.85 

Chandra 

382.85 

(9.37 ±0.56) X 10-13 

(22.04 ±1.32) X 10-13 

(6.33 ± 0.38) X 10^1 

2012 June 10.67 

Chandra 

619.67 

(5.70 ±0.40) X 10-13 

(10.13 ±0.70) X 10-13 

(2.91 ±0.20) X lO-ii 

2012 Oct 7-21 

Swift-XRT 

745.17 ±7.15 

(3.48 ±0.81) X 10-13 

(4.26 ± 1.00) X 10-13 

(1.22 ±0.29) X 10"*! 

2012 Oct 5-Nov 1 

NuSTAR , XMM 750.31 ± 13.33 

(3.64 ±0.17) X 10-13 

(4.00 ±0.18) X 10-13 

(1.15 ±0.06) X lO-ii 

2013 Jan 21-Mar 29 

Swift-XRT 

877.69 ±33.59 

(3.40 ±0.43) X 10-13 

(4.00 ±0.51) X 10-13 

(1.15 ±0.12) X lO-ii 

2013 May 14-Jun 28 

Swift-XRT 

980.01 ± 22.33 

(2.73 ±0.39) X 10-13 

(3.29 ±0.48) X 10-13 

(0.95 ±0.14) X 10-11 

2013 Nov 1.67 

XMM-Newton 

1128.67 

(1.88 ±0.07) X 10-13 

(2.15 ±0.08) X 10-13 

(0.62 ± 0.02) X 10-11 

2013 Dec 11-30 

Swift-XRT 

1178.11 ± 9.57 

(1.50 ±0.49) X 10-13 

(1.72 ±0.56) X 10-13 

(0.49 ±0.16) X 10-11 

2014 Jun 1.25 

Chandra 

1340.25 

(1.53 ±0.13) X 10-13 

(1.82 ±0.15) X 10-13 

(0.52 ±0.04) X loll 

2014 May 11-Jun 27 

Swift-XRT 

1343.01 ± 23.89 

(1.22 ±0.47) X 10-13 

(1.49 ±0.58) X 10-13 

(0.43 ±0.17) X loll 

2014 Nov 30-Dec 24 

Swift-XRT 

1535.60 ± 12.90 

(0.64 ±0.04) X 10-13 

(0.72 ±0.44) X 10-13 

(0.21 ±0.13) X loll 


Note. — Here the fluxes and luminosities are in the 0.2—10 keV range. The fluxes are derived as detailed in The SN explosion date is 

assumed to be 2010 October 1. 


a value the same as that for the 19 keV component, i.e. 
Nh{GSM) = 3.74x 10^^ cm“^, the best temperature goes 
to very low values kT <C 0.1 keV. In the case of powerlaw 
fits, the power law index becomes negative, and the nor¬ 
malization becomes unphysically high. When we ht the 
spectra by fixing the column density for the extra compo¬ 
nent to be the same as that of host galaxy, the reduced 
improved significantly, from 2.71 to 0.94. However, if we 
let Nff^Soit) vary freely, it did not change significantly 
from that of host galaxy column density value, and the 
reduced did not improve. When we kept the column 
density fixed to the host galaxy value, the 0.2-10 keV 
absorbed (unabsorbed) flux of the extra component was 
(1.58±0.67) X 10-13 gj.g g-i ((i.99±o.17) x lO’^^ 
erg cm”^ s”!). If we let the column density vary, the 0.2- 
10 keV absorbed (unabsorbed) flux of the component is 
(1.55±0.14) X 10-13 erg cm-^ 3-1 ((1.93±0.17) x 10-i3 
erg cm-^ s-i). 

If we assume that the column density of the extra 
component is fixed to Nh {Host) in both the 2011 Oc¬ 
tober and 2012 June data, then the component is much 
stronger at the later epoch. Surprisingly this component 
is not present in the 2014 June Chandra data. If we let 
the column density be a free parameter, then the flux 
is roughly constant. We do not understand a process 
that can increase the flux by almost a factor 5 between 
2011 October and 2012 June and then vanish completely. 
Thus we consider the varying column density model to 
be the more realistic model for this extra component. As 
per our fits, the component appeared around 2011 Oc¬ 
tober, then became optically thin (reduced absorption 
consistent with the host galaxy absorption), and finally 
disappeared. The fits to 2011 October and 2012 June 
data are plotted in Fig. [5l 

In the 2014 June Chandra data the SN component is 
fitted with the absorbed apec model and the best fit col¬ 
umn density is Nh{GSM) = (6.82+3;03) x lO^^ cm-^. 
Initially we let the temperature be a free parameter, to 
check if the forward shock cooled down significantly or if 
the reverse shock too has started to contribute towards 


X-ray emission. But the best fit temperature again hit 
the model upper limit of 80 keV, so we fixed it to kT = 19 
keV as in previous datasets. The best fit absorbed apec 
model gives a reduced x^ = 1-15. The extra component 
seen in the 2011 October and 2012 June data is no longer 
detectable. The best fit models are listed in Table IH and 
the best fits for all the Chandra data except for those 
from 2011 October and 2012 June are shown in Fig. HI 
The best fit models for the 2011 October and 2012 June 
data are plotted in FigO 

3.2.3. XMM-Newton data 

For the XMM-Newton observations, we use the best 
fit models of UGC 5189A and the contaminating sources 
from the 2014 June Chandra data (nearest in time to 
XMM-Newton observations) to remove the contamina¬ 
tion in SN flux estimation. We again use an absorbed 
the rmal p lasma model to fit the SN spectra as described 
in 113.2.21 The data are best fit with Nh (GSM) = 
(2.64lg;5g) X 10^^ cm-^ with a reduced = 1-41- We 
plot the spectrum and contour plots in the lower panel 
of Fig. H 

3.2.4. Swift-XRT data 

The Swift-XRT observations are usually closely spaced 
in time but with short (< 10 ks) exposure times. To 
increase the detection significance, we combine several 
data sets to extract the spectra near-simultaneously. 
Our criterion to group the spectra was to get a uni¬ 
form coverage on a logarithmic scale. We, therefore, 
group the observations taken during 2010 November 5- 
20, 2010 November 23-December, 2011 April, 2012 Octo¬ 
ber, 2013 January-March, 2013 May-June, 2013 Decem¬ 
ber, 2014 May-June, and 2014 November-December. The 
spectra were extracted using the online Swi f t-XRT spec¬ 
trum extraction tool (|Evans et al.l 120091: iGoad et "ni 
1200711 . We fitted the apec model with a fixed temper¬ 
ature of 19 keV and C-statistics were used. In 2012 Oc¬ 
tober data (10.7 ks), 2013 December data (14.6 ks), and 
2014 May-June data (13.1 ks), there are only 67, 46 and 
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Figure 7. In the top panel, we plot the 0.2-10 keV light curve 
for SN 2010jl (Blue filled circles). We plot the bolometric light 
curve taken from [Wansson e t al.l <12 0141^ in blue s quares. We also 
plot the X- ray light curve llChandra et ^I2012H) and bolometric 
luminosity JStritzinger et aTTI^O^) of another well studied Type 
Iln supernova, SN 2006jd (orange circles and squares, respectively). 
Unlike SN 2010jl which has a very steep X-ray decay with index of 
—2.12, the SN 2006jd decline is very fiat (index —0.26) in a similar 
time range. In the lower panel we plot the evolutio n of the CSM 
column density for SN 2010jl as explained in Section 13.21 


35 counts, respectively, so we fixed Nh (CSM) to that of 
the best fit values closest in time, obtained from analysis 
of data of other telescopes. The Swift-XRT spectra are 
plotted in Fig. [H 

In the 2010 November Swift-XRT data (Fig. [6|), 
there is a clear indication of an extra feature around 6 
keV, which was also seen in the Chandra data around 
the same epoch. A Gaussian is best fit with an en¬ 
ergy of Acauss = 6.391 q 23 ^he unabsorbed 0.2- 

10 keV fluxes in the Gaussian and the continuum com¬ 
ponents are (3.39 ± 1.28) x 10“^^ erg cm“^ s“^ and 
(26.15 ± 6.67) X 10“^^ erg cm“^ s“^, respectively. In 
2010 Dec data we attempted to add a Gaussian, since it 
was also seen in the 2010 December Chandra data. How¬ 
ever, the data have very few counts thus the fits statistics 
do not change significantly. 

The best fit SN 2010jl models from our analysis are 
listed in Table [H The 0.2-10 keV flux values are given 
in Table m In Fig. 13 we plot the evolution of the 0.2- 
10 keV luminosity as well as the column density. For the 
first 300 days, both quantities evolve slowly, but after 
300 days, one can see a steep power law decline (decay 
index ^ 2) in the luminosity. 

4. RADIO ANALYSIS 

Before analyzing the SN 2010jl data, we examine the 
radio emission in the SN FoV in Very Large Array (VLA) 
archival data taken during 2006 December 1-21 in 1.4 


GHz band. The telescope at the time of these observa¬ 
tions was in G-configuration. The duration for all the 
datasets ranged from 30 minutes to 3 hours (including 
the calibrator). We carried out a joint analysis of the 
data. The image resolution obtained was 16" x 14" and 
the map rms was 191 /rJy (Fig. [5]). The figure shows that 
SN 2010jl lies in a region of extended radio emission. To 
investigate the nature of the radio emission, we looked 
into the VLA NVSS (resolution 45") and FIRST (reso¬ 
lution 5") images of the SN FoV in the 1.4 GHz band. 
In Fig. [SI we overlay the Second Palomar Observatory 
Sky Survey (POSS 11) J-band gray image with the NVSS 
and the FIRST contours. The figure suggests that the 
extended emission is resolved at higher resolutions and 
not likely to contaminate the SN flux. 

The first radio data of SN 2010jl were taken on 2010 
November 6.59 UT in EVLA C-configuration at 8 GHz 
band. The total duration of the observation including 
overheads was 30 minutes. The data quality was good 
and only 7.5% of the data were flagged. The map rms in 
the 8 GHz band was 24 /rJy and the synthesized beam 
size was 2.35" x 2.07". We did not detect SN 2010jl. 
The flux density at the SN 2010jl position was —41 ± 24 
/rJy. We then attempted to observe the SN 2010jl in 
33 GHz band to account for the scenario in which the 
SN was absorbed at lower frequencies. The observations 
were taken on 2010 November 8.47 UT for 3588.7 s. We 
obtained a rms of 45 /iJy and image resolution of 0.67" x 
0.60". We did not detect SN 2010jl in this band either. 
The flux density at 33.56 GHz at the SN 2010jl position 
was 52 ± 45 ^Jy. 

We continued to observe SN 2010jl at regular intervals. 
The first detection of the SN came on 2012 April 18 in the 
22 GHz band, with a flux density of 60.9 ± 17.6 ^Jy. To 
estimate the flux density of SN 2010jl in all the images, 
we fit two Gaussian models, one for the SN component 
and one for the underlying background level to take care 
of any underlying extended emission. Since then we have 
been detecting SN 2010jl in various VLA bands (Table 
Our most secure detections are in 2012 December, 
when the VLA was in the A-configuration. 

To determine the SN 2010jl position, we have used 
2012 December 2 data in 22 GHz band when the VLA 
was in A-configuration. In this data, we obtained a 
resolution of 0.09" x 0.08". The best SN position is 
a = 09''42"*53.®32773±0.00021,,5 = -t09°29'42."13330± 
0.00344 (J2000), which agrees well within 0-15" accu racy 
with the optical position given bv lOfek et al.l ([20141 1. 

The post detection observations of SN 2010jl with the 
VLA G-configuration in 2013 June were contaminated by 
the extended flux due to poorer resolution in this con¬ 
figuration, especially in the 5 and 8 GHz bands. In this 
case, we have used the G-configuration images of 2012 
February before the SN 2010jl detection in the respec¬ 
tive frequencies and subtracted it from the post detection 
images to get the uncontaminated flux of SN 2010jl. 

In Fig. 13 we plot the contour plots of SN 2010jl in the 
5 GHz, 8 GHz and 21 GHz bands for this epoch. The 
SN flux densities obtained in 2013 Jan when the VLA 
was in A—>-D configuration are not very reliable due to 
contamination from the underlying extended emission. 

In Fig. [13 we plot the light curves of the SN (upper 
panel), and in the lower panel plot the spectra at 3 epochs 
when the SN was detected in multiple bands. For the 
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Figure 8. Upper Panel: Pre-explosion SN 2010jl FoV image with the VLA in the 1400 MHz band. The flux units in the upper color bar are mjy. 
SN 2010jl is at a position a. — 09^42^53^.337, (5 — -|-09'^29^42.^^13 (J2000). The immediate SN 2010jl region and UGC 5189A are indistinguishable 
as the image resolution is 16^^ X 14^^. Lower panel: The overlay NVSS (45^^ resolution; left side) and FIRST (5^^ resolution; right side) images 
on the J-band POSS II (Palomer Transient Survey II) image with the same contour levels. Much of the extended emission is absent in the higher 
resolution image. 
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Table 6 

VLA radio observations of SN 2010jl 


Date of 
obsn. (UT) 

Days since 
Explosion^ 

Config. 

Central 

GHz 

BW*’ 

Preq, GHz 

Resolution 

Flux density‘s 

MJy 

rms 

M-Jy 

2010 Nov 06.60 

37.60 

c 

8.46 

0.26 

2.35 

X 

2.07 

< 71.1 

23.7 

2010 Nov 08.48 

39.48 

c 

33.56 

0.26 

0.67 

X 

0.60 

< 133.5 

44.5 

2010 Nov 08.52 

39.52 

c 

22.46 

0.26 

0.93 

X 

0.85 

< 74.4 

24.8 

2010 Nov 14.52 

45.52 

c 

22.46 

0.26 

0.93 

X 

0.86 

< 82.8 

27.6 

2010 Nov 23.43 

54.43 

c 

22.46 

0.26 

1.06 

X 

0.92 

< 70.1 

23.4 

2011 Jan 21.26 

113.26 

CnB 

7.92 

0.13 

3.71 

X 

0.84 

< 76.8 

25.6 

2011 Jan 21.26 

113.26 

CnB 

4.50 

0.13 

6.52 

X 

1.51 

< 116.7 

38.9 

2011 Jan 21.39 

113.39 

CnB 

22.40 

0.26 

0.98 

X 

0.78 

< 107.4 

35.8 

2011 Jan 22.26 

114.26 

CnB 

22.46 

0.26 

1.03 

X 

0.37 

< 79.8 

26.6 

2011 Jan 23.27 

115.27 

CnB 

4.50 

0.13 

5.26 

X 

1.63 

< 76.2 

25.4 

2011 Jan 23.27 

115.27 

CnB 

7.92 

0.13 

2.93 

X 

0.91 

< 55.5 

18.5 

2011 Apr 22.25 

204.25 

B 

22.46 

0.26 

0.41 

X 

0.28 

< 114.0 

38.0 

2011 Jul 07.99 

280.99 

A 

4.50 

0.13 

0.60 

X 

0.42 

< 78.0 

26 

2011 Jul 07.99 

280.99 

A 

7.92 

0.13 

0.27 

X 

0.22 

< 62.1 

20.7 

2012 Jan 24.28 

481.28 

DnC->-C 

8.46 

0.26 

3.35 

X 

2.46 

< 123 

23.4 

2012 Feb 28.20 

516.20 

C 

8.55 

1.15 

2.65 

X 

2.16 

< 32.1 

10.7 

2012 Mar 02.22 

519.22 

C 

5.24 

1.54 

3.89 

X 

3.53 

< 33.6 

11.2 

2012 Apr 15.03 

563.03 

C 

8.68 

1.41 

2.77 

X 

2.25 

< 45.6 

15.2 

2012 Apr 15.05 

563.05 

C 

5.50 

2.05 

3.89 

X 

3.29 

< 32.7 

10.9 

2012 Apr 18.03 

566.03 

C 

21.20 

2.05 

1.13 

X 

0.98 

60.9 ±17.6 

10.1 

2012 Apr 22.10 

570.10 

C 

21.20 

2.05 

1.08 

X 

0.87 

38.3 ±20.7 

13.5 

2012 Aug 11.71 

681.71 

B 

9.00 

2.05 

0.84 

X 

0.71 

76.6 ±20.3 

12.9 

2012 Aug 12.71 

682.71 

B 

5.50 

2.05 

1.26 

X 

1.08 

111.9 ± 17.8 

12.6 

2012 Dec 01.46 

793.46 

A 

5.50 

2.05 

0.34 

X 

0.30 

131.3 ±22.1 

11.1 

2012 Dec 02.40 

794.40 

A 

9.00 

2.05 

0.24 

X 

0.22 

118.8 ±16.8 

9.8 

2012 Dec 02.53 

794.53 

A 

21.20 

2.05 

0.09 

X 

0.08 

115.3 ± 17.1 

9.9 

2013 Jan 18.34 

841.34 

A^D 

21.20 

2.05 

0.17 

X 

0.05 

< 109 

10.7 

2013 Jan 18.38 

841.38 

A-j-D 

9.00 

2.05 

0.32 

X 

0.12 

< 470 

20.9 

2013 Jun 10.96 

984.96 

C 

8.68 

2.05 

2.31 

X 

2.03 

123.0 ±26.8 

9.7 

2013 Jun 11.94 

985.94 

C 

5.50 

2.05 

3.98 

X 

3.23 

91.3 ±34.8 

12.7 

2013 Aug 10.74 

1045.74 

C 

21.20 

2.05 

1.09 

X 

0.90 

< 69.3 

23.1 


^ Assuming 2010 October 1 as the explosion date iStoll et al.ll20i^ 

^ Bandwidth of the observation 

^ Since SN is off the Galactic plane, the errors in the SN flux due to calibration errors will be less than 5%. 
In case of non-detections, the 3-<t flux density limit is quoted. 
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Figure 9. VLA A-configuration 5 GHz, 8 GHz and 21 GHz de¬ 
tections of SN 2010jl, taken on 2012 Dec 01.46, 02.40 and 02.53, 
respectively. The units of the color bars above the maps are fiJy. 


spectra with the JVLA data, we have divided the 2 GHz 
bandwidth into 2 subbands and imaged it independently 
to get the flux densities in the two subbands. 

5. DISCUSSION 

Using the NuSTAR data, we derive a shock tempera- 
ture of 19 keV, w hich is consistent with the analysis of 
lOfek et ^ (|2014ll . This corresponds to a shock velocity 
of ~ 4000 kms“^ at around 750 days. At earlier epochs, 
the shocked gas may be hotter, but we cannot determine 
the exact temperature at any other epoch due to the ab¬ 
sence of hard X-ray observations. This may introduce 
some errors in column density estimate and the SN flux. 
We attempt to quantify this error. In a standard SN- 
CSM interaction model, the temperature of the shock 
varies with time t as (where n is the power-law 

index of the ejec ta density profile), i.e. for a typi¬ 

cal value n = 10 (|Chevalier &: FranssoDl[200^ . Thus over 
the full span of our observations (day ^50 to ^1500), the 
temperature would change by a factor of ^ 2. In our fits, 
we estimate the change in column density and in SN flux 
due to the change in temperature. We find that for a 
20% change in temperature, the column density changes 
by 5 % and the SN flux changes by 1.5 %. If we change 
the temperature of the shock by 100 %, the change in 
column density is less than 15 % and in the SN flux is 
less than 10 %. This analysis shows that even though 
we do not have a handle on the temperature at various 
epochs and are assuming a constant temperature, the er¬ 
rors introduced in other parameters like column density 
and the SN flux do not suffer significantly. 

5.1. Summary of main results 

SN 2010jl is the only Type Iln SN for which a well 
sampled X-ray dataset exists. Thus we are able to trace 
the evolution of column density and light curve all the 
way from 40 to 1500 days (Fig. [7]). 

The X-ray luminosity for SN 2010jl is roughly con¬ 
stant for the first ^ 200 days with a power law index 
of 0.13 ± 0.08. Unfortunately there are no data between 
200 to 400 days, but from day 400 onwards one can see 
a faster decay in the luminosity following a powerlaw in¬ 
dex of —2.12 ± 0.13. For a shock velocity of 4000 km 
s“^ estimated above, the start of the rapid decline in X- 
rays corresponds to a radius of 1.3 x 10^® cm. In Fig. [7l 
we overplot the b olometric luminosity li ghtcurve of SN 
2010jl taken from iFransson et al.l (j2014f l on the X-ray 
light curve of the SN. We note that the bolometric lumi¬ 
nosity also declines rapidly around day 30 0 , cons istent 
with the X-ray lightcurve. IFransson et al.l l|2014ll have 
argued that the steepening in the bolometric luminos¬ 
ity coincides with the Ha shift becoming constant and 
have expla ined it a result o f less efficient radiative ac¬ 
celeration. lOfek et ahl (|2014ll argue that this is the time 
when the SN reached the momentum conserving snow¬ 
plow phase. However, as we will show in Section I5T^ this 
can be explained in our simple model by changes in the 
CSM density profile. 

In Fig. [71 we also plot the X-ray and bolometric lu- 
minosities of another well sampled Type Hn, SN 2006jd 
([Chandra et al.1 l2012bl : iStritzinger et al.l 1201^ . The lu¬ 
minosity decline is much flatter in SN 2006jd than that 
of SN 2010jl. Around the same epoch, the SN 2006jd 


































































Figure 10. Radio light curves of SN 2010jl are shown in the upper panel. The inverted triangles are Scr upper limits. The lower panel 
shows the radio spectrum at 3 epochs. 


X-ray lightcurve declines as while SN 2010jl de¬ 

clines as . The relative flatness of the bolometric 
luminosity of SN 2006jd for a longer duration indicates 
that the CSM interaction powered the light curve for a 
much longer time than SN 2010jl . This indicates that 
the duration of mass ejection in SN 2006jd may have 
been longer in this case than for SN 2010jl, though in 
both cases occurred shortly before the explosion. This 
suggests different nature of progenitors for the two SNe. 

The most interesting result of the paper is the evolu¬ 
tion of the column density with time. At t ^ 40 days, 
the column density is 3000 times higher than the Galac¬ 
tic column density, and declines by a factor of ^ 100 by 
the epoch of our last observation in 2014 Dec. Since the 
higher column density is not associated with the high 
host galaxy extinction, this indicates that the higher col¬ 
umn density is due to the CSM in front of the the shock 
where the dust is evaporated, thus is arising from the 
CSM. 

The presence of bro ad emission lines seen in early SN 
201011 optical spectrajSmith et al.ll2012t Ihdansson et al.l 

12014 iCTek et al.ll2014f) can be explained by electron scat- 

tering (lChugaill2001ll . This requires an electron scatter¬ 
ing optical depth >1 — 3, i.e. a column density greater 
than 3 x 10^"^ cm“^. Comparable values of the column 
density were seen only in the first X-ray observations 
{t < 70 d; Table IH). At later epochs {t > 70 d) we 
only see the X-rays which pass through a column density 


< 10^^ cm“^. These constraints are difficult to reconcile 
with a spherical model of the column density ~ 3 x 10^^ 
cm“^. We therefore should admit that either the width of 
emission lines at the late time is not related to Thomson 
scattering, or the X-rays escape the interaction region 
avoiding the CSM with a high column density. Thus, a 
possible scenario for the CSM in SN 2010jl is th e same 
as th e bipolar geometry in the CSM of 77 -Carinae (|Smi thi 

Hool. 

An intriguing issue here is the presence of an extra 
component in the 2011 October and 2012 June Chan¬ 
dra data, which is not present before or after. This com¬ 
ponent is well ht with a power law index of 1.7, though 
with a varying column density, with much higher column 
density at the 2011 October epoch. The origin of this 
component is not clear. However, the component seems 
to spatially coincide with the the position of the SN and 
occurs when there are changes to the SN spectrum in the 
energy range close to that of the extra component. This 
would seem to suggest that the emission is related to the 
SN. One possibility is that the soft component is a result 
of a cooling shock. In this scena rio a mass loss rate of 
O.lM 0 yr“^ (jFransson et al.l [201411 . wind velocity ^ 100 
kms“^ (jFransson et al.ll2014[ l. and shock temperature of 
19 keV (this paper) corresponds to a cooling time of 86 
days at 10 ^® cm (roughly the radius at one yr for ejecta 
velocity of 4000 kms“^). Thus the forward shock should 
be cooling around this time, as mentioned in our discus- 
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sion of the optical light curve. An adiabatic shock model 
underestimates the low energy X-ray flux. The fact that 
the need for an additional component disappears could 
be a result of the shock becoming adiabatic. 

While the radio observations of SN 2010jl started as 
early as ~ 40 days post explosion, the first radio detec¬ 
tion was around day 566. Unlike its X-ray counterpart, 
the radio luminosity from SN 2010jl is weak for a SN Iln 
(Fig. [Toll . This is unlike SN 2006jd which was an or- 
der of magnitude brigh ter in radio bands than SN 2010jl 
(jChandra et al.l[2012bt) . Like most of the known radio 
Type Iln SNe, SN 2010jl rises at radio wavelengths at late 
times, most likely due to absorption by a high density of 
CSM, or due to internal absorption le.g.. iChandra et ^ 
l2012bri . We have not attempted to fit a detailed model 
to the radio data in view of the small number of de¬ 
tections and small range of time. However, the existing 
data have distinctive features as seen in Fig. [TOl The 
light curves are fairly flat, as are the frequency spectra. 
For the standard models of radio SNe, these properties 
occur when the SN is making the t ransition from op¬ 
tically thick to optically t hin (e.g., iWeiler et al.l 1200^ 
iChevalier fc FranssonI 12?)?)^ . Thus, most likely we have 
detected the radio emission near the peak of the syn¬ 
chrotron emission. 

From Fig. [TUI we estimate that the light curve peak at 
5 GHz occurs on about day 900 at a flux of Fob ~ 0.12 
mJy, or a luminosity of 3.4 x 10^^ ergs s“^ Hz“^. The 
time of maximum is typical of SNe Hn, but the luminos¬ 
ity is lower than most, suggesting that the absorption 
mechanism is not synchrotron self absorption if the radio 
emitting region is expanding at 40 00 —6000 km s~^ as in- 
dicated by X-ray (this pap er and iChaiidra et al.ll2012a[l 
and near IR observat ions (iBo rish et a l.ll2015l) . This can 
be seen in Fig. 1 of iChevalied ( ~09t) . In the case of 
the Type Hn SN 1986J, the radio expansion was mea¬ 
sured by VLBI techniques and an expansion velocity of 
5700 ±1000 kms“^ was found over the period 1999-2008 
([Bietenholz et al.]l2010D . Thus an expansion velocity of 
5000 — 6000 kms“^ is plausible. In this case, the absorp¬ 
tion would have to be something other than synchrotron 
self-absorption, most likely free-free absorption. The fact 
that the 8 GHz radio flux appears around day 700 im¬ 
plies the emission measure along the line of sight at this 
stage {nl£) ^ 8 x 10^® cm“®, where £ is the linear size of 
the absorbing gas region and assuming Tg = 10^ K. On 
day 700 the radius of the shell with an average velocity 
of 4000 km s“^ is of r ^ 2.4 x 10^® cm. Adopting £ ^ r 
we thus come to the rough estimate of the column den¬ 
sity Nh ~ 4 X 10^^a;“^ cm“^, where x is the hydrogen 
ionization fraction. A low hydrogen ionization fraction 
of X ~ 0.1 is needed for the value of Nh to be consis¬ 
tent with the column density on day 700 recovered from 
the X-ray data (3 x 10^^ cm“^). Although we expect 
the CSM to be ionized at early times, recombination is 
possible at later times; however, a detailed study of the 
ionization balance in the CSM of SN 2010jl is beyond the 
scope of this paper. 

Another possibility for explaining the escape of the ra¬ 
dio emission is that it comes from a different region than 
the X-ray emitting region. To have a lower column, the 
radio region would probably come from expansion into a 
lower density part of the CSM and would thus be more 


extended. The X-ray emission from the radio region may 
not be detected because of the low density. Again, details 
are beyond the scope of this paper. 

5.2. Circumstellar interaction modeling 

The detailed evolution of the absorbing column density 
derived from the X-ray data provides us with an unprece¬ 
dented opportunity to examine the CSM around a SN Hn 
in both X-ray and optical bands. The question arises 
of whether the optical light curve powered by the CSM 
interaction is consistent with the observed column den¬ 
sity. Here, we present a simple circumstellar interaction 
model which suggests that freely expanding SN ejecta 
collide with the dense CSM. In a smooth dense CSM the 
interaction zone consists of a forward and a reverse shock 
along with a cool dense shell (CDS) formed in-between. 
We confine ourselves to the interactio n hydrodynamic s 
based on the thin shell approximation (jChevalieilll982t) . 
The equations of motion and mass conservation are inte¬ 
grated for arbitrary density distributions in both ejecta 
and CSM using a 4th order Runge-Kutta scheme. The 
model provides us with the CDS radius {Rg) and velocity 
{vs), the forward shock speed Vfs ^ Vg, the reverse shock 
speed Vrg = {R/t — Vg), and the kinetic luminosity re¬ 
leased in the shock Lj = 2TTR‘lpjVj, where pj is preshock 
density and Vj is the shock velocity (j = rs, fs). 

Generally, the conversion of the kinetic energy into 
the radiation in the interaction SNe is affected by com¬ 
plicated hydrodynamic a nd thermal pro cesses including 
the thin shell instability (|Vishniaclll994[) . the Rayleigh- 
Taylor instability of the decelerating CDS, the CSM 
dumpiness, mixing, and energy exchange between cold 
and hot components via radiation and thermal conduc¬ 
tivity. This makes the computation of the radiation out¬ 
put of the CS interaction quite a formidable task. We use 
a simple approach in which the X-ray luminosity of both 
shocks is equal to the total kinetic luminosity times the 
radiation efficiency rj = t/(t + tc), where C is the cooling 
time of the shocked CSM at the age t. To calculate the 
shock cooling time we assume a constant postshock den¬ 
sity four times the upstream density po, while the shock 
temperature is calculated in the strong shock limit. We 
find that the reverse shock is always fully radiative in our 
models. In the case of SN 2010jl, the optical luminosity 
dominates the observed X-ray luminosity by a factor of 
ten. In our model, the bolometric luminosity, which is 
primarily in the optical, is equal to the total radiation 
luminosity. 

The model assumes that the optical radiation gener¬ 
ated in the interaction zone instantly escapes the CSM 
which means that the diffusion time (tdif ~ (r/c)r, 
where c is the speed o f light and and r is the shell 
radius, iChevalien I198H1 is smaller than the expansion 
time, or the CSM optical depth t < cjv (where v is 
the shell speed). In our model, on day 3 the column 
density of the wind is 5 x 10^® cm“^ and the Thomson 
optical depth at this stage is therefore 30. This value 
is comparable to c/v ~ 30 assuming v = 10^ kms“^. 
Therefore, after about day 3 the diffusive trapping of 
photons in the wind can be ignored. The SN ejecta 
density distribution is approximated by the analytical 
expression p^ oc 1/[1 ± (x/uo)"] which reflects an in¬ 
ner plateau {v < vg) and an outer power law density 


























ACTION IN SN 2010jl 


19 



Figure 11. Top Panel. The bolometric light curve obtained from 
our modeling detailed in {dashed line) pl otted over the observa - 

tional data {blue squares) which are taken from f^ansson et al.l 11201411 . 
The red solid line is the shell velocity obtained from our modeling (in 
units of 10^ km s~^) compared to the observational data {filled circles) 
taken at various epochs (this paper on day 750 and IBorish et ^12013 , 
at earlier epochs). Lower Panel: the CSM column density outside the 
forward shock obtained from our modeling in JEl versus the data re¬ 
covered from X-ray observations {squares). Our simple model detailed 
in El is able to reproduce the observational quantities quite well. 

drop with 7 < n < 10. The initial ejecta boundary 
velocity Vb is assumed to be 3 x 10^ km s“^; our re¬ 
sults are not sensitive to the boundary velocity in the 
range of (1 — 3) x 10^ km s“^ at epochs t > 5 d. The 
CSM density distribution is set by a broken power law 
p oc with m = mi = 2 in the range of r < ri and 

m = m 2 > 2 for r > ri. The value mi = 2 was cho¬ 
sen for small radii because it is an approximate fit and 
is plausible for a wind; the outer value of m 2 is a fitting 
parameter. The ejecta diagnostics based on the CSM 
interaction cannot constrain the ejecta mass and energy 
uniquely because the same density versus velocity dis¬ 
tribution in the ejecta outer layers can be produced by 
a different combination of mass and energy. Yet the ob¬ 
served interaction luminosity and the final velocity of the 
decelerated shell (u/) constrain the energy of the outer 
ejecta with the velocity v > Vf. For the adopted density 
power law index n the obvious relations vq oc 
p{v) oc po{vo/v)'^, and pQ oc M/v^ result in the energy- 
mass scaling E oc or FI oc assuming 

n = 8. In turn, this scaling, when combined with the re¬ 
quirement that the velocity Vf should be larger than the 
ejecta density turnover velocity vq, provides us with the 
lowest plausible values of E and M. As a standard model 
we adopt an ejecta mass of 8 Mq which barely satisfies 
the condition Vf > vq. The corresponding kinetic energy 
of ejecta is then 2.1 x 10®^ erg. 

The optimal model is found by the y^-minimization 
in the parameter space of ri, m 2 , and the wind den¬ 


sity parameter wi = Anr^p in the region r < ri. We 
found the best-fitting values wi = 4.1 x 10^^ g cm~^, 
ri = 1.29 X 10^® cm, and m 2 = 4.45. The parameter 
errors determined via variation do not exceed 1.5 %. 
These errors are formal and too optimistic, since we ig¬ 
nore uncertainties in the distance and the extinction, and 
the systematic error related to the model assumptions. 
The model (Fig. [TT]) reprodu ces the SN 2010jl bo lomet- 
ric light curve (obtained from iFransson et ^1201411 quite 
well and produces reasonble values for the shell veloc¬ 
ity, consistent with the observa tional estimates (o btained 
from this paper on day 750 and IBorish et al.l|2015L for the 
earlier two epochs). The evolution of the column density 
obtained from X-ray data is also described by our model 
at t > 200 d as well, except for the early epoch t < 100 
d when the model requires a somewhat larger column 
density. In the model, the CSM density power law index 
breaks at ri = 1.3 x 10^® cm from m = 2 to m = 4.45, 
which suggests that the bulk, 2.6 Mq, of the total CSM 
mass, 3.9 Mq, lies within the radius ri. In order to de¬ 
scribe the column density at late epoch t > 600 d, the 
power law index should become steeper (m ~ 7.5) for 
r > r 2 = 3.2 x 10^® cm. A lower density region out- 
si de the close-in CSM is consistent with what is deduced 
bv iFransson et all (j2014[ l in their analysis for the narrow 
lines. 

The agreement between the optical model CSM col¬ 
umn density and that inferred from X-ray data suggests 
that the CSM density recovered in the model is realistic. 
The external radius of the CSM envelope, 1.3 x 10^® cm, 
combined with the wind velocity of 100 km s“^ implies 
that the SN event was preceded by an episode of vigor¬ 
ous mass loss starting 40 yr prior to the explosion. The 
mass loss rate at this stage was ^ 0.06 Mfn yr~^ , whic h 
is similar to the value obtained in iFransson et all (1201411 . 

Although the overall picture of the X-ray generation by 
CSM interaction with absorption in the CSM is generally 
convincing, there is an issue with the X-ray luminosity. 
The point is that the unabsorbed X-ray luminosity is sig¬ 
nificantly (a factor of 10) lower then the optical. But at 
late times, t > 100 d, the forward shock dominates the 
luminosity so one would expect that more than half of 
the total X-ray radiation escapes the interaction zone, 
at odds with the observations. The disparity suggests 
there is some mechanism that converts the kinetic lu¬ 
minosity into the optical radiation avoiding significant 
hard X-ray {hi/ > 1 keV) emission. The soft XUV radia¬ 
tion then could be absorbed by the CSM resulting in the 
X-ray deficit. Given the Thomson optical depth tt < 1, 
hvjmt? <C 1, and AkTf,/mc^ <C 1, neither Compton scat¬ 
tering in the CSM nor Compton cooling of hot electrons 
in the forward shock are able to provide the requir ed de¬ 
gree of X-ray softening (|Chevalier fc Irwiiil 1201211 . An 
alternative scenario is conceivable which connects the 
X-ray deficit to CSM dumpiness. Although radiative 
properties of the shocked CSM cannot be reliably quan¬ 
tified, the predominance of soft radiation is a likely out¬ 
come in this case. Indeed, if the bulk of the CSM were 
in clumps the luminosity of the shocked intercloud gas 
would be weak. The cool matter mixed with hot gas in 
this scenario becomes a dominant source of the radiation 
that falls into the soft XUV band. Thermal conductivity 
might allow heat flow from the hot shocked intercloud gas 
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into t he mixed cool gas. Alternatively, iFransson et al.l 
(|2014f) have explained this deficit to be due to the pres¬ 
ence of an anisotropic CSM, because of which most of the 
X-rays are absorbed and converted into optical photons. 
We note that in the case of SN 2006jd, the bolometric 
luminosity is around an order of magnitude larger than 
the X-ray luminosity (Fig. [7]). 

5.3. Comparison with Other Results 

In our X-ray spectral fits of UGC 5189A and other 
nearby sources, we allow the column density to vary 
freely and also use a metallicity of 0.3 (iStoll et al.ll201lll 
for the excess column density (over Galactic). We de- 
rive a host colum n density of 4 x 10^^ cm“^, whereas, 
lOfek et all (|2014ll have fit all the nearby X-ray sources 
with a fixed Galactic column density (3 x 10^° cm“^) 
and a power law spectrum with photon index F = 1.375. 
Part of the discrepancy can be accounted for the fact 
they have used solar metallicity as opposed to 0.3 so¬ 
lar metallicity for the host Galaxy used by us. This 
will lower their equivalent Hydrogen col umn density in 
the X SPEG fits by a factor of three. IFransson et al.l 
(IMl . using a fit to the Lyman-a damping wings, find 
Nniiiost) = (1.05 ± 0.3) x 10^° cm“^ for the host 
galaxy, while the corresponding value is A^jj(Galactic) = 
(1.75 ± 0.25) X 10^° cm“^ from the Milky Way. Their 
Galactic column densi ty is around a factor o f 2 low er than 
the one derived from iDickev &: LockmanI l)1990ll . How¬ 
ever, the GMRT has observed the SN 2010jl host galaxy 
UGC 5189 in 21 cm radio bands (|Chengalur et al.ll20l4) 
in 2013 November-December, and they derive the HI col¬ 
umn density to be 2.4 x 10^^ cm“^, consistent with our 
best fit values. We caution here that the column densi¬ 
ties obtained in our fits have large uncertainties, and our 
estimates for the host galaxy may be treated as an upper 
limit on the host column density. 

Althou gh our bes t fit t emperature agrees with that 
found bv lofek et al.l (j2014l l. our column density is some¬ 
what smaller {Nh{CSM) = (6.67l^ g 4 )x 10^^ cm“^) than 
that quoted bv lOfek et al.l (I2014D which is ^ 10^^ cm“^; 
the difference is more significant considering they have 
used solar metallicity. To test the robustness of our best 
fit column density we have attempted to fit only the 
XMM-Newton spectrum as XMM-Newton data are espe¬ 
cially sensitive to the column density due to absorption 
being dominant in the 0.2-10 keV energy range. We ob¬ 
tain a column density of Af_y(CSM) = (6.67l4 gg) x 10^^ 
cm“^ {x^ = 1.07), which is consistent with our joint 
fit. The contour plot of the column density for XMM- 
Newton data shows that it is well constrained (Figure 
^ _ 

lOfek et ^ ()2014f l dispute the values of high column 
density at early times obtained bv iChandra et all (j2012ari 
by claiming that they have used many parameters. How¬ 
ever, in addition to the fits to X-ray data, the addi¬ 
tional evidence of a high column density comes from the 
Fe K-a line seen in 2010 December 7-8 Chandra data 
and 2010 November Swift-XRT data, which suggests 
Nh{GSM) = 2 X 10^^ cm“^ (assuming ZjZQ = 0.3 and 
EW = 0.2 keV), consistent with the absorption column 
d ensity obtai n ed in our fits. 

lOfek et ^ (120141) fit the early Chandra data with 
blackbody models, considering the medium to be opti¬ 


cally thick to X-rays. However, a blackbody fit corre¬ 
sponds to an extremely small emitting area. For ex¬ 
ample, at an early epoch, for an X-ray luminosity of 
8.5 X 10"^^ erg s“^ and a temperature of 3.4 keV, the 
blackbody emitting radius is only 2 x 10^ cm, which is 
physically not plausible, considering the expected large 
area of the shock front. T hus we disfavor the blackbody 
model. lOfek et ^ (1201411 also fit their models using a 
powerlaw. In our models, we have tried to fit the data 
using powerlaw models. While a powerlaw does give ac¬ 
ceptable fits, the photon index is very flat (F < 0.5), 
which is physically implausible. 

The steepening in the X-ray and bolometric lumi- 
nosity ligh t curv e around day 300 is quite significant. 
lOfek et ^ (|2014ll explain the steepening by arguing that 
the shock reached the fast cooling snow-plow phase. 
However, in our model we can explain this by intro¬ 
ducing a steepening in the density profile. Our mass 
loss es timates are a factor of 10 smaller than lOfek et al.l 
(IMll . but are consistent with IFransson et al.l (|2014l) . 
This discrepancy can be partly accounted by the fact 
that lOfek et ^ (j2014f l have assumed a wind velocity of 
300 kms“^, whereas we assumed i t to be 100 kms“^, 
adopted from IFransson et all (|2014ll . 

6 . CONCLUSIONS 

Here we have reported the most complete X-ray and 
radio observations of a luminous Type Hn supernova SN 
2010jl. SN 2010jl is the only Type Hn SN which has been 
well sampled in both radio and X-ray bands since early 
on. 

Using publically available NuSTAR data, we deter¬ 
mine a temperature for the shocked gas that is consis¬ 
tent with the value obtained from lOfek et ^ (j2014l) . 
The 6.4 keV Fe-Ka lin e seen in the first Chandra epoch 
(jGhandra et al.l l2012a| ) is also present in the Swift- 
XRT data taken around the same time, confirming that 
this line is real. 

While the radio emission is weak in SN 2010jl , the 
X-ray luminosity is one of the highest for a Type Hn SN. 
This provides us a unique opportunity to trace the evo¬ 
lution of the circumstellar column density. The circum- 
stellar column densities at various epochs are 10-1000 
times higher than that of the host galaxy. This evo¬ 
lution is satisfactorily reproduced in a model in which 
freely expanding SN ejecta collide with the dense and 
smooth GSM, with the forward shock luminosity being 
a fraction of the kinetic luminosity depending upon the 
radiation efficiency. 

The X-ray light curve evolution is quite fiat for the 
first 200 days. However, after ^ 400 days, the light curve 
shows a steep decline. There is a similar steepening of the 
optical luminosity, from which we infer a steepening of 
the GSM density power law index from the standard 
profile. In contrast, the Type Hn SN 2006jd light curve 
at a similar epoch is much flatter, indicating differences 
in the mass loss leading up to the SN. The case of SN 
2010jl is consistent with rapid mass loss beginning a short 
time before explosion. 

The observed radio emission for SN 2010jl is very weak 
and does not clearly evolve as in standard models. The 
radio spectra and their evolution suggest that the emis¬ 
sion is close to its peak at an age of ^ 10^ days. The 
implication is synchrotron self-absorption was probably 
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not a factor in the rise to maximum and that another 
process, likely free-free absorption, dominated. 
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